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Abstract 

The aim of this work is to study the holographic dual to the gauge theory with a nonzero gluon 
condensate. We check for consistency the holographic way of describing the condensate and calculate 
the expectation value of a small Wilson loop in the presence of the gluon condensate, thus obtaining 
the relevant coefficient in the operator product expansion of the small loop in different holographic 
models. We also study the effect of the condensate on the Gross-Ooguri phase transition in the 
correlator of two circular Wilson loops in parallel and concentric configurations. In the numerical 
study of the concentric case, we find that the phase transition changes its order when the size of the 
loops is of order of the gluon condensate. We report this change of the phase transition order to be a 
new effect in Wilson loop correlators. 
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I. INTRODUCTION 



Holographic models provide a powerful machinery for the nonperturbative description of 
gauge theories in the strong coupling regime. Formulated first as a method of description of 



the strongly coupled supersymmetric conformal gauge theory N=4 SYM jjfl, the AdS/Cr- 
correspondence was then used in a phenomenological approach to study the real QCD [4|. Its 
deformations breaking the supersymmetry and conformal symmetry were proposed as well in an 
effort to consistently reduce the AdS/CFT correspondence to the description of the a QCD-like 



theory 



One of the main manifestations of the broken conformal symmetry in QCD is the nonzero 
gluon condensate, the nonperturbative vacuum expectation value of the scalar gluon operator 

(«,tr(G 2 )), (1) 

introduced in the framework of QCD sum rules in 8|. This value plays a major role in the 
QCD sum rules method and has a variety of phenomenological applications. For instance, one 
can show, that the p and other meson masses are proportional to (a s tr(G 2 )) 1//4 , and the energy 
of QCD vacuum equals e = — ^(^tr(G 2 )). Therefore, there are strong reasons to believe 
that in the real QCD the vacuum expectation value under consideration does not vanish. The 
phenomenological value of the gluon condensate which is usually used is of order of 

(a s ti{G 2 )) ~ (200 MeV) 4 . (2) 

In the framework of the QCD sum rules the gluon condensate is a free parameter, which can 
be adjusted in order to get better phenomenological results. Later, in the 1980-ies, in order to 
determine whether the condensate really exists, it was measured in the fundamental theory via 
lattice calculations. For this purpose one can study the expectation value of a small Wilson 
loop operator 

W(C) = ±r{tr Pexp[<tigA,dx»]). (3) 



As was shown in a sufficiently small Wilson loop can be presented in the form of an operator 
product expansion, which looks like (omitting the input of the perturbation theory) 

,y( C ) = 1 -i<^!»^ + ..„ (4) 



where S is the area of the Wilson contour C. This relation was a basis for the lattice study of 
the gluon condensate, which has led to a reasonable value, while the procedure of subtracting 



the perturbative contribution into the Wilson loop is not well defined |10l. 

In the gauge/gravity duality the gluon condensate can be calculated in two at first sight 
unrelated ways. First of all, as the gluon condensate is a vacuum expectation value of a certain 
operator, it can be extracted as a coefficient in front of the normalizable mode of its dual field 
in the supergravity In fact there is a general way to show that the field dual to the tr(G 2 ) 
operator is the dilaton jl2 |. That means that to describe a field theory with a nonzero gluon 
condensate we need to consider the supergravity background with a nontrivial normalizable 
mode of the dilaton field 

= 00 + 04^, (5) 

where z is a distance to the boundary of the holographic background. 

On the other hand, one can study the area of the minimal surface, spanned on the circular 
contour of a small radius on the boundary, which is related to the expectation value of the 
field-theoretic Wilson loop on this contour, and get the value of the gluon condensate in the 



spirit of lattice calculations, as it was done in 



13j. In this approach the gluon condensate is 



extracted from the OPE of the Wilson loop operator. In fact, the relation (j3J) is valid only in the 
small coupling limit, while the holographic calculation implies a strong coupling. That means 
that we need the version of the Wilson loop OPE at a strong coupling, which was obtained 
in [l^] in the AdS/CFT ap proach. Most of the results from 14] were later checked in gauge 
theory calculations 15|, lis], so we treat them as a strong coupling supersymmetric equivalent 
of 01]). As a result of this calculation we discover the literal coincidence of the two ways of 
describing the gluon condensate in holographic models. 

Looking for the nonsupersymmetric version of the Wilson loop OPE we find that its form is 
rigidly determined by the dimension of the stack of the Dp-branes, which form the background of 
the model. Thus it has the same form in all models, based on D3-branes, being insensitive to the 
IR behavior of the model, and it is different in the Sakai-Sugimoto model p, H); surprisingly 
coinciding with the small coupling result fl4]). We note here, that the applicability of our 
approach to QCD is rather subtle, because the OPE is valid at a small enough size of the loop 
(namely r <C (a s tr(G 2 )) -1 / 4 ), which implies a large enough energy scale. On the other hand, 



the supergravity description is valid only when the coupling constant is large, which in the 
case of QCD implies a small energy scale. Recall thought, that in QCD the coupling already 
becomes large at the scale of lGeV, which means that there is a window in the energy scale 
where our approach is applicable. 

In the second part of the paper we study the effect of the the gluon condensate on the 
Gross-Ooguri phase transition in the correlator of two circular Wilson loops. In a configuration 
with two small parallel loops of the same radius we evaluate the shift of the phase transition 
point due to the condensate. In a coplanar concentric configuration we rely on the numerical 
analysis and find an interesting change of the phase transition order, when radii of the loops 
become close to (a s tr(G 2 ))~ 1//4 . 

The work is organized as follows: in the first section we briefly review a simple holographic 
model that describes a nonzero gluon condensate and calculate the correction to the expectation 
value of the Wilson loop due to the presence of the condensate. In the second section we 
discuss the form of the OPE of a small Wilson loop in a general p-brane background and in 

n 

the D4 Witten background |7j in particular. Sections three and four are devoted to the study 
of the Gross-Ooguri transition in the presence of the condensate for parallel and concentric 
configurations, respectively. In the Appendix we provide the normalization procedure of the 



dilaton field in the model and the calculation of the gluon condensate, which is similar to 18 



but is included here for the sake of consistency and is rewritten in a slightly different notation. 



II. ONE LOOP, SUPERSYMMETRIC CASE 



We start by considering the holographic model based on the background of a stack of N c 
D3-branes. Its most famous feature is that in the throat limit this geometry looks like the 
AdS^ x S 5 space and exhibits a superconformal symmetry, the same as the N=4 SYM theory 
The supersymmetry can be partially broken by D7-branes, representing the fundamental 



quarks 



19] , the conforma 



scale parameters 



(3, Q, Eo 



symmetry of AdS$ can be violated by the introduction of various 



2l|, but is restored in the vicinity of the boundary. As QCD is 



asymptotically free in the ultraviolet region (dual to the vicinity of the AdS boundary), we 
wish for the holographic dual of QCD to have this kind of boundary behavior. To study the 



gluon condensate in the model we need to take the supergravity solution with a nontrivial 
dilaton background (jSJ). A simple background, involving a nonzero gluon condensate is that of 
a D-instanton smeared on the D3-branes 2l| . This background in the string frame has a metric 



ds 2 m \str = - T y/hZ'i{dx tl dx^ + dz 2 + z 2 dQ 2 5 ) (6) 



z 2 



and a dilaton 



g 9 h- u = 1 + Y W 



where A is related to the 't Hooft coupling constant of the dual field theory. The curvature radius 
of the AdS is related to the string coupling g s (the dilaton asymptotic value) as L 4 = A7rg s N c l 4 , 
where l s = \fa! is a string length scale. We see, that in the large N c limit this space is slightly 
curved and the supergravity description is valid. The dimensional parameter q violates the 
conformal symmetry which is restored on the boundary z — > 0. 

The string coupling constant g s is related to the coupling constant of the dual Yang-Mills 
theory as follows \12\. 



9ym = Aix 9s- (8) 

To proceed with our study of the gluon condensate (tr(G 2 )) we need to fix the relation 
between the dimensional parameter q and (a s tr(G 2 )) . This procedure has been performed in 
1^ | . but we repeat it in the Appendix in a different notation, used in this paper, for the sake 
of consistency. The result is flA.8j) : 

i = ^<a,t r( g)) 
' \ s/2\ N c 1 ' 

Now we can find the analogue of (J4]) in this model. In order to compute the Wilson loop 
via holography one needs to study the minimal area of the string worldsheet spanned on the 



Wilson contour at the boundary of the bulk space [22] 



(W{C)) = e - Area{c \ (10) 

The area of the string worldsheet that we are interested in is described by the Nambu-Goto 
action in the string frame 



Sng= 2^P 1 



where g is the induced two-dimensional metric on the worldsheet. In the string frame metric 
of the D3-model the action (1111) for the circular Wilson contour of the radius R takes the 
form (we expand the metric in 4 , because 04 ~ ^-qcd anc ^ we are considering the worldsheet 
in the region of the bulk space z ~ R <C ^qq D ) 

r 

2vr 



.s 







SNc= 2^ dr ^ rv ^+ T ( 1 + ^ + ) • ( i2 ) 



We assume here the parametrization of the worldsheet by the radius and angle variables and 
perform the integration over the latter due to the symmetry of the problem. At this stage 
we end up with a one-dimensional problem of finding a solution to the equation of motion, 

n 

following from (1121) . This solution at 04 = has a simple form [15( 



Zo(r) = VR 2 - r 2 , (13) 

and after the subtraction of a linear divergency gives a constant contribution to the area of the 
minimal surface. We can calculate the correction due to the gluon condensate to the Wilson 
loop expectation value using the perturbation theory in 04-z 4 , and then plug this correction into 
the action (I12jl in order to get in the linear order in 04^ 4 (the superscripts denote the order in 



J C°(z + z x ) + C\z ) = J jC°(z ) + C\z ) 



S_£P_ 

5dz Zl 



(14) 

boundary 



The last term in the integral vanishes as it is the equation of motion, and the boundary term 
must vanish as well, because by definition the boundary value of the correction is 0. However, 
it does not in reality, because the expression SC S qI°^ diverges on the boundary when calculated 
on the zero-order solution, thus giving a finite input to the result. One can subtract this 
divergence by introducing some boundary counterterms, but we will use a somewhat more 
elegant approach. We shall simply rewrite everything in a different coordinate system, where 
this divergence do not exist at all. 

The form of the zero-order solution f|T3|) gives us a hint on what coordinate system to choose. 
We will use the spherical coordinates, where the solution looks simply like a constant. We define 
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them as 

* = /sin(0), (15) 
r = /cos(0), (16) 

so that the zero-order solution is simply fa(9) = R. The Nambu-Goto action ( !T2|) in these 
coordinates is 

tt/2 

s «° = I / d9 7»F ^WT7(¥ (i + ■ (it) 



One can check that the boundary term in f|14jl vanishes because &C & q^ = 0. Thus we end up 
with a simple result coming from the full action, evaluated on the zero-order solution. 

tt/2 

L 2 f cos(#) 04-R 4 . ,_. 2 /10 . 

o 

The first term, as expected, gives the usual linear divergency, proportional to the perimeter of 
the contour, while the second term defines the condensate part of the Wilson loop OPE: 

A/HA li2 ^ B 4 ^ (aMG 2 )) p4 

where we have used the relation ([9]) and the definition of L. Note that we get the right sign 
of the correction and right A^-dependence as compared with (j4j). Although the dependence 
on the coupling constant is not the same as in (J4]), this is not very surprising, because this 
result is computed in the strong coupling limit. If we compare our calculation of the Wilson 
loop in the presence of the gluon condensate with the computation of the OPE coefficients of 

n n 

a Wilson loop from [14|, we find a literal coincidence. In [14[ the OPE is computed via the 
exchange of a dilaton between two spherical loops. It is presented as the integration of the 
dilaton vertex operator over two surfaces, multiplied by the propagator of the dilaton field. 
The vertex operator is nothing but the correction to the action due to the dilaton that we 
have denoted C 1 , and the propagator gives a factor z 4 , which is present in the profile of the 
dilaton in our model. It turns out that by calculating the first order action on the zeroth order 
solution, we calculate the same integral as in Thus what we get does coincide with the 
strong coupling version of the Wilson loop OPE, and we conclude that the two ways of defining 
the gluon condensate in holography are in fact identical. 
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III. ONE LOOP, NONSUPERSYMMETRIC CASE 



To describe QCD holographically one needs to introduce a gravitational background, where 
supersymmetry and conformal symmetry are broken. As we have already mentioned, in the 
models based on D3-brane background [5|, |20J, the OPE of the Wilson loop will look like (fT9|) . 
3ut these models are not the only case. For instance, the Sakai-Sugimoto AdS/QCD model 
17] is based on the D4-brane background, developed in Q. For the sake of generality, we will 



consider the case of N r p-branes here. For an arbitrary p this supergravity solution in the string 



frame looks like 



23N25| see 12[ for a review) 



dsi 



str 



h 2r] ab dx a dx b + (du 2 + u 2 dVt\ 



(20) 



with a dilaton 



where in the throat limit 



g s h(u) P 4 3 : 



h{u) 



7-p 



L 7 ~v = d p g s N c l 7 ~ p , 



d p = (4tt) Vr 



7 — p 



(21) 

(22) 

(23) 
(24) 



and indices a,b denote the (p+l)-dimensional space parallel to p-brane. We can make the 



coordinate transformation 



26] 



u 



P 



2 

"5-p 



L 5 ~p z 5-p 



(25) 



which brings the metric to a conformally flat form similar to ([6]) 



str 



7-p 



5 — p L\ 5 -p 



r] ab dx a dx b + dz 2 



P 



z 2 dn 2 8 _ p 



(26) 



Given this expression we can easily compute the induced metric on the string worldsheet with 
a circular boundary and hence the Nambu-Goto action (fTi]) . 



S 



NG\Dp 



dr 



7-P 

5 — p L \ -'-p 



rV z 2 + 1 



(27) 



We note, that the string coupling constant g s as well as N c enter this expression only via the 
parameter L. Making this dependence explicit, we find 

Sng\d p ~ (g s N c )^ ~ A ^ (28) 

We see now, that the dependence of the Nambu-Goto action on y/\ is rigidly related to the 
dimension of the underlying D-brane worldvolume. Let's consider the case of p=4. The string 
frame metric of this setting in the conformal coordinates ( 1251) is 
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ds 2 str = (J^J ^ u dx»dx u + f(z)dr 2 + j^dz 2 + Az 2 dnl\ . 



(29) 



Here the r dimension is compactified with the period 5r = ^-Zo, the black hole warp factor is 
f(z) = 1 — ^ and /i, v stand for usual 4D space indices. The dilaton is 

3 



e* = 9s ( #V (30) 



2z) 

The parameter L is related to the string coupling constant as 

L 3 = 7tg s N c l 3 s (31) 
The relation between g s and the field theory coupling constant can be established similarly to 



the D3 case 



27] 



671 i s 



28 



for 



Let us note here, that the metric (|29|) is not conformal from the very beginning (see 
a review), as the rescaling of 4D coordinates can't be absorbed now in the rescaling of z. It is 
not asymptotically AdS and its curvature behaves as -, making the classical gravity approach 
inapplicable near the boundary. Therefore the quantum gauge theory dual to this geometry is 
not asymptotically free. This is a significant obstacle for us, because now our method of the 
dilaton operator normalization ( 1A.4|) is useless. Indeed, in an asymptotically non-free theory 
the consideration of the leading divergence of a perturbation theory ( IA.5j) is not legitimate even 
in the limit of asymptotically large momenta. Thus, if we consider a perturbation of the dilaton 
profile by a dimensional quantity $4 (as in the preceding section), we cannot unambiguously 
fix its relation to the gluon condensate. 

9 



Nevertheless, let us assume, that in QCD a simple relation between the dilaton profile and 
gluon condensate, similar to (Q, holds at least parametrically. In this case, we can proceed in 
our calculations. 

The Nambu-Goto action in the D4-brane background fl29l) is (from now on we omit all 
numerical factors, as we are interested on parametrical dependence only), see ([28]) . 

R 

SngIda ~ {g a N c )l a J^y/l + f{z)z* (33) 



Note that the scale of compactification z , enters now the equation of motion for z(r) via the 
warp factor f(z). We do not know the solution to the classical equation of motion, but we can 
analyze it from the point of view of dimensional analysis. The integral in ( 13"3"|) has a dimension 
(-1) in the units of length. There are two scales, R and zq, that can contribute to this dimension, 
but in the limit R — > the integral must remain finite in order for it to agree with the leading 
term in fll]). Therefore the value of the integral computed on the classical solution must have 
the form —g(R), where g(R) has no negative powers of R. If we introduce now the parameter 
q/X of dimension (-4) in the dilaton profile, then, similarly to (1191) . it will enter this function 
as g(R) = const + jR 4 + 0(q 2 R 8 ). At the end of the day we can find the parametric structure 
of the gluon condensate correction to the small Wilson loop 

6{W) ~ (g s N c ) ^(0 4J R 4 + ...)- Mr A (G2)) S 2 , (34) 

where we have used the expressions for g s (1321) and 4 (jHJ). Interestingly, it is exactly the 
behavior that has been presented in (T4J). 

IV. GROSS-OOGURI TRANSITION, PARALLEL LOOPS 

When studying the impact of the dilaton on the Wilson loop vacuum expectation values it 
is quite natural to expand the consideration to Wilson loop correlators. Particularly it is very 
instructive to investigate the behavior of a correlator of two parallel Wilson loops of equal radii 
R that are separated by a distance / in a transversal direction, which we shall do in this section, 
and of concentric coplanar Wilson loops of different radii, which we are going to study in the 
next section. 
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As it has been proposed by Maldacena |22[, the Wilson loop correlator may be calculated in 
string theory and is related to the surface area of a string worldsheet stretched over two Wilson 
contours. Generally speaking, two types of the worldsheet solution are classically admissible - 
a connected and a disconnected one, and we have to choose the one that has the minimal action 



(and thus the surface area) of the two. As it has been first discussed in [29|, when distance 
between the contours increases and reaches a certain critical value I — l c , a disconnected 
worldsheet solution becomes more preferable than the connected one: 

^conn.ij'c) ^disc- 

Thus the connected Wilson loop correlator undergoes a second order phase transition - the 
Gross-Ooguri phase transition - and vanishes. (Actually, the disconnected surfaces can still 
exchange propagating gravitational modes, the account of which transforms the phase transition 
into a crossover. However, these effects are suppressed by N c and will not be considered in our 
paper.) At an even greater value / = the connected solution becomes unstable and ceases to 
exist. An analogous transition is possible in the case of two concentric Wilson loops, when it 
is determined by the ratio of their radii. 

These effects in the pure AdS^ x S 5 targe t sp ace in the two possible configurations as 



described above have been analyzed in detail in 30( and 3l| respectively. In the case of parallel 



loops that will be under our consideration in this section the critical value equals l c = 0.91-R 
and the connected solution becomes unstable at 1* = 1.04R. The same transition has been 



studied in a nonconformal background numerically in [32|. In this section we shall determine 
how the particulars of the Gross-Ooguri transition in the case of parallel loops are modified by 
the presence of a gluon condensate. 

A. Equations of motion and boundary conditions 

Let us consider a pair of Wilson loops in a Euclidean R A space, parametrized as 

d ={{x 1 ) 2 + (x 2 ) 2 = R 2 , x 3 = -1/2, x° = 0} 
C 2 ={(x x ) 2 + (x 2 ) 2 = R 2 , x 3 = 1/2, x° = 0}. 
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Let us also denote x 1 = r cos tp, x 2 = r sin tp, x 3 = x, x° = t. The correlator of the Wilson loops 
is related to the Nambu-Goto action corresponding to the minimal surface stretched between 
the two Wilson loops in the AdS^ x S 5 space with a dilaton 22[: 



(W(C 1 )W(C 2 ))=ex V (-S 



NG, 



The metric tensor of the AdS$ part has the form 



ds 2 = (1 + 4 ^ 4 ) 1/2 ^{dz 2 + dt 2 + dx 2 + dr 2 + r 2 dtp 2 ). 
Due to obvious symmetries, the minimal surface will be parametrized as 

t(r,o-) = 0, (f(r,a) = a, x(t, a) = x(r), r(r, a) = r(r), z(r, cr) = z(r). 



This defines the action 



5 



NG 



^Jdr^{z' 2 + r' 2 + x' 2 ) l/ \l + ^) 



4x1/2 



(35) 



where a prime denotes a derivative with respect to r. A minimal surface is one that satisfies 
the equations of motion for the action (135jl . The form of the action implies that the momentum 
k along the x-coordinate is conserved. 

We shall find the minimal surface solution perturbatively in 04. 



B. Zero-order solution in different coordinates 



We will first discuss the zero-order solution (zq(t), J"o(t), xo(t)), i.e. the solution without the 
dilaton. In the analysis below we will follow somewhat closely 30]. Let us choose first r = x. 
The zero-order equations of motion yield: 



( r o + z o + x o)xx — 0- 



(36) 



Let us recall the boundary conditions: 



z{-l/2) = z(l/2) = 0, r(-Z/2) = r(l/2) = 0, 
2^(0) = 0, r'(0) = 0. 



(37) 
(38) 
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These boundary conditions are the same both for the full and for the zero-order solutions. The 
latter (138j) reflects the fact that the solution is symmetric with respect to the x = hyperplane. 
As a result, we obtain from (I36| [37J [38]) that 

zl + rl + x 2 = R 2 + l 2 /4 = a 2 . (39) 

The coordinate x is denoted here as Xq for the sake of universality with respect to different 
coordinate systems. 

Now we shall choose different coordinates, namely, = r, and 

z = /(0)sin0, r = /(0)cos(0), sc = a;(0). (40) 

Due to the symmetries of the solution, it is clear that the angular coordinate runs from to 
some value 0o as x varies from —1/2 to 0, the region < x < 1/2 corresponds to running back 
from O to 0. The boundary is located at = and we shall determine 8 later. In the new 
coordinates the action of a connected worldsheet may be rewritten as: 

do 

S conn . = 2 V\ [ d6 -^L (f + f + .,/») V* (i + 04 f S in 4 0) 1/2 , (41) 
J j sin 
o 

where a prime denotes now a derivative with respect to 0. 
Zero-order equations of motion for the action (1411) are: 



/ o sin 2 X ' 



U/o 2 + /o 2 + ^)~ 1/2 = ^ (42) 



COS %/o(/o 2 + /o 2 + -o 2 )- 1/2 



/ o sin^0 



COS0 , 2 /2 /2 \l/2 COS0 / 2 /2 #o\ -1/2 ,. n% 

2— 2^(/o+/o 2 + ^) +— 2^(/o+/o 2 + *o) • (43) 



/q sin sin 

Eq. (|42l) allows us to express x through f : 



»_ k*f*sm*9{f* + f») 
cos 2 0-P/ o 2 sin 4 ' 

The boundary conditions are transformed into: 

x(0) = -1/2, x(6 ) = 0, (45) 
x'(0o) = oo, (46) 
/(0) = i2, /' (0 O ) = 0. (47) 

13 



The last boundary conditions in (14T)} H7|) may be obtained from the analysis of the behavior 
of d6/dx, df /dx in cylindrical coordinates. Also note that the zero-order relation ( 1391) implies 
that /o + x o = const = a 2 . This means that we can replace the boundary conditions (|47|) with 
the equivalent conditions that are more suitable for numerical calculations: 

/(to) = a, f'(8 ) = 0. (48) 

Using Eq. ( 144]) and boundary condition ( 1461) we obtain the expression for 9q through k: 



# = arccos = [ka\. (49) 

\ 2ka I 

Substituting x' from 044p into (1451) we get a single nonlinear second-order differential equation 
for fo(6) with two boundary conditions ( 14H|) . We shall solve it numerically, with the solution 
fo(8;l/R,ka) depending on two dimensionless parameters: known l/R and for the moment 
unknown ka. If we substitute the numerical solution fo(8; l/R, ka) into Eq. (1441) . we may now 
integrate it, and taking into account boundary conditions (1451) . we obtain: 

e [ka] 







1/2 = x (9 ) - x (0) = k dd x l l^M-±I^ [l/R, ka]. (50) 

cos 9 — k z jQ sin 



This allows us to determine the yet unknown ka for each value of l/R and find the solution 

Ue-i/R). 

C. The on-shell action 

We can define the zero- and first-order corrections to the action (J4TJ as follows: 

S° conn . = 2V\[d6 (f + + 1/2 , (51) 

J /osm 6» 
o 



0o 

SLn. = ^04 | d0 / 3 cos 61 sin 2 61 (/ 2 + /^ 2 + x' 2 ) 1/2 





X CO r cO 

T h\P) + 73 x i(^J 



Sf J y ' 5x 

where (fo,x ) and (fi,xi) are the zero- and first-order solutions respectively. 

An important property of the solution f (9;l/R) is that its derivative f' Q {6) ~ 0{6 2 ) when 
6 — > 0. This means that the boundary terms in (152j) 
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-^f^ = const, = 0, ^2^ = k. (53) 

The first order solutions have to satisfy the following boundary conditions: 

zi(0) = x x %) = A(0) = 0. (54) 

Combining (I53[ 154)) we conclude that the boundary terms in (|52|) turn out to be zero. Hence, 

substituting x' from (jUJ) into f )52l) we get: 

fo / \ 1/2 







SL. = / rf^ /o 3 cos 9 sin 2 (9 I /° 2 + JLg ) ■ (56) 

q \ J cos 2 6 / 

Let us now make a comment on the choice of coordinates. In the papers 3CH32j the calcula- 
tions are carried out with r = x. In our case staying in those coordinates would lead to nonzero 
boundary terms in ( |52|) (see discussion in Section II). Moreover, those terms would contain a 
singularity proportional to the dilaton, which cannot be cancelled out by subtracting a term 
proportional to the contour perimeter. In order to make sense of the singular expressions in 
those coordinates we would have to add boundary counterterms in the spirit of holographic 
renormalization. The choice r = 9 allows us to deal with finite and regular first-order correc- 
tions straightaway. 

D. The Gross Ooguri phase transition: numerical results 

The Gross-Ooguri phase transition occurs when the Wilson contours are separated by such 
a great distance l c that the disconnected solution becomes preferable to the connected one: 

Sconn.iJ'c) £>disa (57) 

where Sdisc. is the action on the disconnected solution ( TTgj) : 

tt/2 



SL, = 2^X [d9 ^ 
J sin 9 



o 
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The factor 2 here stands for the two connected parts of the worldsheet surface. If we separate 
now the actions and the distances into zero- and first-order quantities, we obtain: 

Sconn.(lco) = ^disci (59) 
i _ ^disc. ~ ^conn. (jeo) (f\n\ 

Cl ~ dS° con Jdl(l c0 ) ■ m 
Note that I c q is determined from Eq. ( 1591) as a distance at which non-singular quantity 
Sconn.(lc) — iSdisc. equals zero. Indeed, using expressions in (1551 1551) it may be rewritten as: 
QO _ cO _ _ 2 /T 1-singp 

conn. disc. " • fl 

Sill fo 

/— /" 1 /' 2 cos 2 # + k 2 f£sm 4 6 , N 

2v /A / rf» — _ ^— ^ -— . (61) 

* 8111 /o (cos 2 (£) - P/ 2 sin 4 0) 1/2 f(/ 2 + /^ 2 ) 1/2 + /o (l - F/ 2 ^|) 

We have already noted that fb(0) ~ 0(6* 2 ), 6* — > 0. Hence, the integrand in Eq. ( 1ST]) is finite at 
8 = 0, and, although it is mildly divergent at 8 = 8q, the integral is finite. This implies that the 
derivative — jp^) = c °""^ — —(Ida) is regular as well. Substituting f Q (0;l/R) into 
the actions (!55ll6Tl) for various ratios l/Rwe can numerically determine that S® onn (I) — S® isc = 
with a margin of error 10 -6 when I = l c0 = 0.9047 R and the derivative in (160]) equals: 

f)q0 

" (lco) = (1-20 ± 0.01)VXR-\ (62) 



conn. 



dl 

as well as find the first-order correction to the action on a connected solution: 

Slonnilcv) =0.26 VA^. (63) 
Combining (HB1 EDI EH E3J we obtain: 

l cl = (6.05 ± 0.06) x 10~ 2 <p 4 R 5 . (64) 

Thus we discover that the presence of the gluon condensate drives the point of the Gross-Ooguri 
transition to larger distances between the loops. An interesting question is how the presence 
of the dilaton affects the point of the onset of the instability Z*. However, it is impossible to 
determine it in a perturbative analysis. In the next section we shall investigate the concentric 
loop configuration numerically. 
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V. GROSS-OOGURI TRANSITION, CONCENTRIC LOOPS 



To study the correlator of two concentric coplanar Wilson loops in the presence of the gluon 
condensate we cannot rely on the perturbation theory analysis, because in this case we cannot 
simply take a small enough radius of the loops to ensure the inequality 04-R 4 <C 1. Thus we will 
treat this case numerically. In this section we study the correlator of two concentric circular 
Wilson loops with radii R\ and R2, lying in the plane x = 0. While we treat this problem 
only classically, other interesting stringy effects were studied in a similar setting in 33j and the 
semiclassical approximation was inspected in a recent paper 34j. Generally speaking, there 



are two possible forms of the minimal surface, ending on this boundary: two disconnected 



half-spheres and one connected surface, similar to a half-torus. As it was pointed out in [311 ]. 
the connected surface exists only when the distance between loop contours is not very large. 
In our numerical study we will discover that this is indeed the case. 

Considering the form of the boundary conditions (two concentric circles on the boundary of 
AdS) and the discussed above divergencies in the (z, r)-coordinates, it will be very convenient 
to adopt for our study the toroidal coordinates (r, a), defined as 

cosh(r) — cos(cr) 

z{r,a) = — — — , (66) 

cosn(r) — cos(cr) 

where a is a parameter of the coordinate system. The coordinate surfaces of constant r = tq 
are the tori with an axial circle of the radius acoth(ro) and the radius of the section circle 
. , a , n , . The coordinate a runs from to 2ir, being a kind of an angular variable in the section 

smh(T u ) ' 

of the tori. We can choose the parameter a in such a way that one of these tori satisfies our 
boundary conditions: 

a = y/RxR 2 , (67) 

cosh(ro) = ?±±*±. (68) 
ri\ — K2 

Hence our task is to find the minimal surface, described by the function r(er) with boundary 
conditions 

r(0) = t(tt) = r . (69) 
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The Nambu-Goto action ffTT]) in toroidal coordinates looks like 

da y/l + T'((T) 



S 



NG 




sinh (r(cr)) 



sin (a) 



sin 2 (ex 



(70) 



v cosh(r) — cos(cr) 

Note that as we are not treating 04a 4 as a small parameter, we do not expand the square 
root factor in the metric We note here that before we have introduced the dimensional 
parameter 04 the theory was conformal, and the only way the dimensional parameter a enters 
the action is in the combination 4 a 4 . Thus we have 2 dimensionless parameters in the problem: 
To fl6Tl) and 



(71) 



We solve the equation of motion, following from ( 1701) . with the boundary conditions (16911 nu- 



merically using the shooting method, obtaining the solutions at different (see Fig, 1(a)). To 
calculate the area of the corresponding surface, we regularize the action by subtracting the sur- 
face area of a torus with the same boundary: r(a) = To- One can check, that the action on such 
a torus contains only the divergent term proportional to the perimeter. Hence by subtracting 
it we remove this divergency from the computed result. 
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FIG. 1. The minimal surfaces with boundary on two circles (a) and one circle (b) calculated at 
= 0, 5, 10, 40, respectively from top to the bottom. 



At = we compare our results with the ones obtained in 31( and find the agreement 
in the case when R\ = 2 R 2 , which is computed there. We also see that at a certain small 
t < t*, corresponding to large R1—R2 ffTT]) . the connected solution does not exist. At = we 
find that r* = 1.405, giving Ri = 2.72 R 2 . The Gross-Ooguri transition takes place when the 



regularized area of the connected surface becomes equal to (—2) - the area of the disconnected 
surface. In the case <fi = we obtain the transition at r c = 1.5357, or Ri = 2.40 i?2- 

To study the transition at <fi ^ 0, we need to solve a nonperturbative equation of motion 
for the disconnected surface as well. As it is just a pair of deformed half-spheres, we use the 
spherical coordinates like in (ITT]) and for different obtain numerical solutions, depicted on 



Fig, 1(b) 
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FIG. 2. The area of the connected (blue) and disconnected (red) surface, depending on tq\ a) at eft = 1 
the Gross-Ooguri transition is smooth, b) at <fi = 6 the Gross-Ooguri transition has a jump. The green 
dashed line represents the minimal area. 



Depending on the value of <fi we can look at two interesting points in tq: the r* where the 
connected solution ceases to exist, and the r c , where the area of the connected solution becomes 
larger than the area of the disconnected one. Obviously, r c cannot be less then r*, because in 
that region there is no connected surface whose area can be measured. The case when r c > 



is depicted in Fig. 2(a) But as we saw in the preceding section, the account of the gluon 



condensate drives the point of the phase transition to a larger separation of the loops (164)) and 
thus to smaller r c , and at a large enough condensate it can reach the value r*. Hence there 
can be a situation when the area of the connected surface remains less then the area of the 
disconnected all the way up to the point, where the connected surface ceases to exist, as it is 



depicted in Fig. 2(b) In this particular case the phase transition happens at r*, and then the 
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change of the form of the minimal surface is accompanied by a jump in the value of the minimal 
surface itself. Thus the Gross-Ooguri transition can change its order depending on the value of 
the gluon condensate. We find, that this is indeed the case and determine at what particular 
value of <p this change occurs. In Fig. |V] we plot r* and r c , determined at different values of 
0. There is a point cr , where these curves intersect. It is located at cr ~ 4.3. Recalling 
our definition of dimensionless parameters flTTj) we find that the characteristic size of the loops 
when the jump occurs is related to the gluon condensate as 

(a s tiG 2 y l/ \ (72) 

Interestingly enough, we find that the coefficient is of order of unity, so we can state that this 
change can indeed be relevant to the usual hadronic physics at the energy scales less than 1 
GeV. 



Null T 

X 

1.5 - 




(a) (b) 



FIG. 3. The dependence of r* (blue) and r c (red) on eft. a) the difference is not visible in this scale, b) 
enlarged dashed region of the first plot, the point, when the change of order occures is visible. 



VI. CONCLUSION 



In this paper we have studied the gluon condensate and its relation to the Wilson loops 
in holographic models. We have calculated the expectation value of a small circular Wilson 
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loop in a supersymmetric holographic setting and found the agreement with the OPE of the 
Wilson loop, calculated in jl^]. Although in the nonsupersymmetric case of the holographic 
model, based on the compactified D4 Witten background 7J, we have been unable to find exact 
results due to the nonconformality of the theory in the UV, nevertheless, we have obtained the 
parametric structure of the OPE coefficient under consideration. It turned out to be the same 
as in the weak coupling calculation {9]. 

Studying the effect of a nonzero condensate on the Gross-Ooguri phase transition in the 
correlator of two Wilson loops, we have discovered that it drives the point of the transition to a 
larger separation between the loops, which is quite natural, since the presence of the condensate 
makes the area of the disconnected surface larger. In the course of a numerical study of the 
correlator of two concentric loops we have found that at values of the condensate comparable 
to the order of the loop radii, the Gross-Ooguri transition changes its order, being accompanied 
by a jump not only in the derivative of the action, but in the area of the minimal surface itself. 
This is a new effect in the Wilson loop correlator behavior and it is worth further investigations. 
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Appendix: Normalization of the dilaton 



Here we rewrite the procedure of the dilaton field normalization and the gluon condensate 
definition in the D3 model, used in 



in the notation of the present paper. 
Using the symmetry arguments we can state that the supergravitational dilaton field is 
dual to the gauge theory operator O^, which is proportional to scalar gluon operator: = 
c^a s ti{G 2 ). Hence, according to the standard AdS/CFT recipe [2|, yj the boundary value of the 
dilaton field is treated as a source of and its normalizable mode as a vacuum expectation 
value (0^). We can compute the two-point function (O^O^) to fix the coefficient c^. In order 
to do this we need to calculate the classical action of the dilaton fluctuation and take the 
second variation with respect to its boundary value. The kinetic term of the dilaton field has 
a canonical form in the Einstein frame metric, which is defined as = yf~g~ s e~^G In the 
Einstein frame the part of the supergravity action that we are interested in, is 

Se = f d XQ x^G~ E (-1^0^ (A.l) 

c 

4(2^) 2 

In the last line we assumed, that there is no dynamics along the S^-sphere, and used the 
definition of L given above. The solution to the equation of motion near the boundary is 

<f>(z, Q) = cf> 9^K 2 (Qz), 0(0, Q) = 0o, (A.3) 

where K 2 is the McDonald function of the second kind. Substituting this solution back into the 
action and taking the second variation we get the expected two-point function. In the leading 
order of the large-Q expansion it is exactly the conformal result 

N 2 1 

<W> = ^gQ 4 m(Q 2 e), (A.4) 

where e is the AdS space cutoff, which is related to the UV cutoff in the quantum field the- 
ory. This expression can be compared with the leading order in the large-Q expansion of the 



/V 2 f 1 

J d A xdz — ((9 2 0) 2 + (^0) 2 ) (A.2) 



35) 



correlator of scalar gluon operators, obtained in the QCD sum rule approach 

N 2 - 1 

(tr(G 2 )tr(G 2 )) = — Q 4 ln(Q 2 e 2 ), (A.5) 
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and this comparison in the large N c limit allows us to fix the operator, dual to the dilaton 

1 



a 



4V2 



(A.6) 



Similarly, we can compute the vacuum expectation value of 0$. As in the considered background 
the classical profile of the dilaton has a normalizable branch, <f> 



jz 4 , we find a nonzero result 



4(2tt 2 



[z,Q) 



d z <p(z, Q) 



N 



«,(2tt 2 



-q 



(A.7) 



Thus, we find the desired expression for the gluon condensate in this model 

(a s tr(G 2 )) = 4v/2a s O = N c ^q 



(A.8) 



Note that the final relation doesn't include a s , and since (a s tr(G 2 )) ~ N c , q does not depend 
on N c . Therefore all the dependencies on N c and a s are reproduced by holography. 
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